An exact solution for the Casimir force in a spherically symmetric medium 
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We calculated the force of the quantum vacuum, the Casimir force, in a spherically symmetric 
medium. Maxwell's fish eye, surrounded by a perfect mirror and derived an exact analytic solution. 
Our solution questions the idea that the Casimir force of a spherical mirror is repulsive — we found 
an attractive force that diverges at the mirror. 

PACS numbers: 03.70.+k, 77.84.Lf 



O 

1—5 



Oh 
I 



> 

cn 

q 

o 



> 

X 



Casimir suggested an intriguing model that could ex- 
plain the stability of charged particles and the value of 
the finestructure constant ^. The argument goes as 
follows: Imagine the particle as an electrically charged 
hollow sphere. Two forces are acting upon it: the elec- 
trostatic repulsion and the force of the quantum vacuum, 
the Casimir force — presumed to be attractive [2HS] • The 
stress a of the quantum vacuum on a spherical shell of ra- 
dius a must be given by a dimcnsionless constant times 
hc/a^ on purely dimensional grounds — the quantum 
stress is an energy density proportional to h, and hc/a^ 
carries indeed the units of an energy density. Now, the 
electrostatic energy of the sphere is proportional to the 
square e^ of its charge and is also inversely proportional 
to a"* [9]. Therefore, an attractive Casimir force balances 
the electrostatic repulsion regardless of how small a is, 
provided e^ /{he) assumes a certain value given by the 
strength of the Casimir force. This strength depends on 
the internal structure of the particle — the fact that it 
is a spherical shell — but not on its size, which could 
be imperceptibly small. Casimir's model, however crude, 
could simultaneously explain the finestructure constant 
e^ /[he) and the stability of charged elementary particles! 

All one needs to do is calculate the Casimir force on a 
spherical shell, but such calculations are notoriously dif- 
ficult. After a marathon struggle with special functions, 
Boyer succeeded in numerically computing the force for 
an infinitely conducting, infinitely thin shell and found a 
surprising result PUI that shattered Casimir's idea: the 
vacuum force is repulsive and so cannot possibly balance 
the electrostatic repulsion. Boyer's heroic calculation 
was confirmed in a sophisticated and elegant paper by 
Milton, DeRaad and Schwinger ^lj and by others [T^ . 
The spherical shell has become the archetype for a shape 
that causes Casimir repulsion [T3]. However, doubts have 
been lingering about whether the repulsive force of the 
shell may be an artefact of the simple model used [l4] . 
for the following reason: the bare stress of the quantum 
vacuum is always infinite and this infinity is removed 
by regularization procedures [2HH]- The most plausible 
regularization involves considering the relative stress be- 
tween or inside macroscopic bodies. But an infinitely thin 
sphere does not represent an extended macroscopic body, 
nor multiple bodies. Suppose the physically relevant vac- 




FIG. 1: The Casimir force on a spherical shell (left) is re- 
pulsive ^10!, or is it? We assumed the shell to be filled with a 
medium (right) and found an attractive force. The shades of 
grey indicate the profile of the medium (plotted in Fig. 2). 

uum stress of an extended spherical shell tends to infinity 
in the limit when the shell becomes infinitely thin and in- 
finitely conducting. In this case the regularization would 
remove this physically significant infinity, producing a fi- 
nite result that may very well have the wrong sign. Our 
paper supports the contention that the Casimir repulsion 
of the spherical shell could be an artefact of regulariza- 
tion. 
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FIG. 2: Index profile n{r) (grey curve) of the medium inside 
the shell and the resulting vacuum stress a{r) (black curve, 
in units of hc/a^). As r — >■ a the stress a tends to — oo. 

Consider a minor modification of Casimir's model 
(Fig. 1). Imagine that the spherical shell (though still 
infinitely conducting and infinitely thin) is no longer hol- 
low, but filled with a medium of gradually varying electric 
permittivity e and magnetic permeability ^. In this way 



we have extended the shell to a macroscopic body where 
the Casimir stress gradually builds up. For preserving 
the spherical symmetry we assume that e and /i depend 
only on the distance r from the center of the sphere. 
We expect zero Casimir force in the center, because in 
a spherically symmetric medium the center does not dis- 
tinguish any direction for a force vector to point to, and 
so the force must be zero. The Casimir stress tensor a 
will be radially symmetric and may change with increas- 
ing r. Considering how a varies we obtain a physically 
well-defined Casimir-force density V • tr. In particular, 
we assume a toy model for e and /^ where a turns out to 
have an exact solution. Our model is (Fig. 2) 
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where tii is a constant. Equation M is valid for r < a, 
and at r = a we place a perfect spherical mirror that sets 
the transversal components of the electric field strength 
to zero. Equation (fTl) describes Maxwell's fish eye [T5] 
used in perfect imaging [l6l - fT8| . We have chosen this 
model because, in our calculation of the Casimir force, we 
take advantage of the mathematical fact that Maxwell's 
fish eye implements the geometry of a simple curved 
space: it represents the 3-dimensional surface of the 
4-dimensional hypersphere in stereographic projection 
(Fig. 3) [19] . As in other applications of transformation 
optics [ini HO] such a geometrical interpretation of elec- 
tromagnetic media can give us guidance for non-trivial 
design problems or, as in our case, for calculations that 
would otherwise be complicated [2T]. We have found a 
simple, exact expression for the vacuum-stress tensor: 
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where 1 is the unity matrix. The stress is isotropic, nega- 
tive and falls monotonically, so the Casimir-force density 
V • (T is always attractive in our model (Fig. 2). Close 
to the mirror the stress and the force tend to infinity. 
Our model thus suggests that the Casimir force at a 
perfect spherical mirror is indeed infinite. An imper- 
fect mirror, on the other hand, may lead to a finite and 
possibly attractive vacuum force, which offers new hope 
for Casimir's fascinating explanation of the finestructure 
constant and the stability of elementary charged parti- 
cles [1]. Of course, in a more realistic theory the particle 
should not be regarded as being a classical object inter- 
acting with the quantum vacuum, as in Casimir's case [T] , 
but rather as a self-consistent quantum structure. Specu- 
lation aside, we found a non-trivial exact solution for the 
Casimir force. Analytic solutions for Casimir forces are 
extremely rare — the attraction between two plates with 
infinite e [2^, the repulsion between plates with infinite 
£ and n |23] and the attractive force on a homogeneous 
spherical ball with infinite £ [21] have been solved; we 




FIG. 3: The medium (plotted in Fig. 2) represents the ge- 
ometry of the hypersphere (shown here as a sphere) in stereo- 
graphic projection (Hnes). The figure illustrates the reflection 
(points) at the mirror (black circle) on the hypersphere and 
in physical space. 

believe we have discovered the first exact solution in a 
gradually varying medium. 

We use Lifshitz theory [31 [S] for our calculation, be- 
cause this is the best physically motivated and tested 
theory of Casimir forces [7]. Lifshitz theory relates the 
vacuum stress to the electromagnetic Green function (as 
in Schwinger's source theory [21] )• The physical stress a 
of the quantum vacuum is expressed as 

a = lim [CT(r,ro) -cro(^,'ro)] , 

1 



CT(r,ro) = T(r,ro) - -TrT(r,ro) 



(3) 



h?Gs{r, Tq] in) dn , 



where the r are the correlation functions of the fields in 
the vacuum state between the points r and Tq at equal 
times; the r are finite for r y^ Tq. The stress is regular- 
izered by subtracting a bare vacuum stress from cr(r, Tq) 
in the limit Tq — >■ r. Our calculation turns out to be 
independent of the actual regularizer ctq a-s long as cto de- 
pends only on local properties of the medium. The total 
correlation function r consists of the sum of the electric 
and the magnetic field correlation functions given by [8] 

hce{r) 
Tel = 

TT 

he [°° ^- 

Tmag = 7 — T / V X G^ (r , To ; Ik) X Vo dn (4) 
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where the arrow indicates differentiation from the right; 
Gs denotes the symmetrized electromagnetic Green func- 
tion (G -|- G"^)/2 for purely imaginary wavenumbers Ik 
{i.e. for imaginary frequencies ick). The Green function 
G describes the electric field at the spectator point r 
generated by a point source at Vq pointing in all possible 
spatial directions. The Green function is thus a bi-tensor 
that obeys the wave equation 



1. 
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£K''G = 6{r-ro)t. (5) 

In our case [Eq. (ITj)] the medium is impedance-matched. 



e(r) = ^{r) = n{r) . 



(6) 



We show that for impedance matching the electric cor- 
relation function equals the magnetic one. For this, we 
represent G as 



G = G 



5{r - ro) 1 



(7) 



and obtain from the wave equation [Eq. ([5|] 

Vx-VxG + nK2G = , , , \ ., (8 

n n[r)n[ro)K^ 

where we expressed the double curl of the delta-function 
term [Eq. (l7|] in terms of derivatives with respect to r 
and Tq. Notice that the magnetic Green function, defined 
as 
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obeys the same wave equation [Eq. ([8|]. Consequently, 
Gmag agrees with G apart from a delta-function term, 
but such a term does not matter in the correlation func- 
tions [Eq. Q] where we regard r ^ Tq before we take the 
limit ro — > r. We conclude that t^i = T,„a„ and obtain 
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(10) 



which sets the scene for calculating the Casimir stress 
a from the electromagnetic Green function G in 
impedance-matched media. 

The Green function for Maxwell's fish eye [Eq. (fTl)] has 
been obtained already in connection with perfect imaging 
[T7] . We state the results we need here. For keeping our 
expressions uncluttered we set 



a — \ , ni = \ 



(11) 



in our calculation and then obtain the general result 
[Eq. &\ by scaling arguments. Suppose the material 
[Eq. (nj] extends to infinity (without the mirror) . In this 
case [18) : 
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(13) 
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In the geometrical picture behind Maxwell's fish eye ^\ 
(Fig. 3) the electromagnetic wave propagates on the sur- 
face of the hypersphere from source ro to spectator r 
in stereographic coordinates with distance arctanr', and 
D denotes the Green function of a conformally coupled 
scalar field [TB] ■ The effect of the mirror is described by 



an adaptation of the method of images [5] on the hyper- 
sphere (Fig. 3). There the mirror lies on the equator (a 2- 
dimensional surface for the 4-dimensional hypersphere). 
We subtract from Go the electromagnetic wave generated 
by the image source on the hypersphere (Fig. 3). This 
field is the mirror image of the original field. In stereo- 
graphic projection [TH], the reflection at the equator cor- 
responds to the transformation r ^ r^^. To obtain the 
reflected wave Gf,, we thus perform the coordinate trans- 



formation r = r{r') with r' 



and then replace r' 



by r. Note that we also need to transform the field com- 
ponents of Go , which is done by the Jacobian ^25: 
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such that 



G^ = PGo(r-i). (16) 

In this way we obtain for the Green function 

G = Go(r)-PGo(r-i). (17) 

One verifies that the transversal components of G vanish 
at r = 1, as they should at a perfectly reflecting electric 
mirror. Note that also for the reflected wave the mag- 
netic Green function equals the electric one (up to an 
unimportant delta-function term) even if the mirror is 
not made by an impedance-matched material. To prove 
this, consider Gmag deflned by Eq. ([9| for the trans- 
formed Green function Gq that describes the reflected 
wave and rename r as r' (we recall that Gq is the re- 
sult of a coordinate transformation) . Then we make use 
of two geometrical facts. First, n{r)~^'V x Go defines a 
one-form with respect to the effective geometry with line 
element ndl [H]. Second, this one- form is invariant un- 



der the transformation r' — 



because the line element 



ndl is invariant. Therefore we can read n(r')~^V' x Gq 
as the coordinate-transformed n{r)^^W x Go- Hence we 
can also read the magnetic Green function of Gq (with 
r renamed as r') as the coordinate-transformed Gmag of 
Go. For Go the entire medium is impedanced-matched, 
and in such a case we have already established that the 
magnetic Green function agrees with the electric one up 
to a delta-function term. Consequently, the same must 
be true for the transformed Green function Gq that de- 
scribes the reflection at the mirror. 

Now we are ready to calculate the Casimir stress from 
the electric Green function according to Eqs. ^ and ( 10 1. 
The Green function Go of the infinitely extended fish-eye 
medium corresponds to the Green function on the entire 
surface of the hypersphere, which is a uniform space. It 
can only produce a uniform vacuum stress ctq that does 
not contribute to the Casimir-force density V • ct. As we 
are interested in contributions to a that do generate a 
force we take the uniform ctq as our regularizer. In this 



way, we are independent of actual regularization proce- 
dures that attempt to explain why the physical vacuum 
stress is finite. As an additional bonus, we only need to 
focus on the reflected part of the radiation field, similar 
to the Lifshitz theory [S] for the vacuum stress between 
conducting plates [5] or in other piece-wise uniform pla- 
nar materials 0]. We thus consider only —PGa{r^^) in 



the total Green function [Eq. (17)]. As the Casimir stress 
in a spherically symmetric medium must be spherically 
symmetric, we calculate —PGo{r~^) only in x direction, 
i.e. we put y = z = and yo = zq = in the definitions 
[Eq. ( [T2p5l )] and evaluate -PGo(r-i) at xq == x-\ We 
obtain, after some straightforward algebra. 
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with the matrix elements 
dD 



di = 2 



dr' 



__dD_ , d^D 
'""d^"'' d^ 



where, in the x direction, r = x and 



1 /l 



(18) 



(19) 



(20) 



Equations (18 [20) enter the formula for the correlation 
function [Eq. ( 10 )] in place of G^ where Gg is integrated 
over all positive-imaginary wavenumbers Ik. It is wise 
to perform the k integration before the r' differentia- 



tions [Eq. (19)]. We obtain for the scalar Green function 
[Eq. (Ill])], using integral 2.4.4.1 of Ref. [26] . 
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which gives 
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167r2 (r'r)'^ 



(21) 



(22) 



All three eigenvalues of r — tq are identical in the x di- 
rection and, by virtue of spherical symmetry, they must 
be identical in all directions. Equation ( [22| ) is thus valid 
everywhere in the medium, and from Eqs. ([3]) and (20) 
we get our result [Eq. (pi)] for a = 1 and ni — 1. Rein- 
stating units for a and r produces Eq. ([2]) for general a. 
For obtaining our result for general ni we notice that ni 
of Eq. (fl]) appears in the wave equation [Eq. ([5])] as uik 
and a prefactor of rii of the source term. Consequently, 



we only need to replace G(iK) by niG{iniK) in Eq. (10) 



where n also carries the prefactor ni, take niK as a new 
integration variable and obtain Eq. ([2]) in full generality. 
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